Uninfected random walkers in one dimension 
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We consider a system of unbiased diffusing walkers (A0 < > QA) in one dimension with random 

initial conditions. We investigate numerically the relation between the fraction of walkers U(t) 
which have never encountered another walker up to time t, calling such walkers "uninfected" and 
the fraction of sites P(t) which have never been visited by a diffusing particle. We extend our study 
to include the A + B — > diffusion-limited reaction in one-dimension, with equal initial densities of 
A and B particles distributed homogeneously at t = 0. We find U(t) ~ [P(t)] 7 , with 7 ~ 1.39, in 
both models, though there is evidence that a smaller value of 7 is required for t — * 00. 



I. INTRODUCTION 

Random walks model a host of phenomena and find 
application in virtually all areas of physics [|l]-[5] . In this 
paper we consider a one-dimensional (ID) system of non- 
interacting, unbiased diffusing walkers A% < — > 0^4, with 
random initial conditions. Recently, we have studied the 
persistence properties of this model [[|. The persistence 
probability P(t) has been widely studied in recent years 
|i] fl9| |. In the present context P(t) is defined as the frac- 
tion of sites which have never been visited by a diffusing 
particle up to time t. Such sites are termed "persistent 
sites". We showed || that, in the case of the A% < — ► %A 
processes, 



P(t) = cxp[-(16p 2 Dt/n) 1 ^} 



(1) 



where po is the particle density and D is the diffusion con- 
stant. The fraction of unvisited sites in this model can 
equally be thought of as the surviving fraction of immo- 
bile particles in the presence of randomly diffusing traps. 
At t = 0, consider each site to contain a static particle, A, 
so that pa{0) = 1. Then for t > 0, let particles of another 
type, B, randomly distributed on the lattice with density 
Po, diffuse though the system with the B particles acting 
as perfect traps for the A particles so that A + B — > B. 
This model is just the well-studied "scavenger" model of 
Redncr and Kang [^-^2| . The fraction of unvisited sites 
is then just pA(t)- 

The focus of this paper is to explore the relation be- 
tween P(t) and the fraction of walkers U(t) that have 
never encountered another walker up to time t. To fa- 
cilitate our discussion, we introduce the following defini- 
tion. We consider all walkers to be uninfected at time 
t = 0. The walkers are, however, considered to be mutu- 
ally infectious, so that for t > any contact between 
walkers leads to mutual contamination. We therefore 
define an uninfected walker to be one which has never 
encountered (occupied the same site as) another walker. 
We can then refer to U(t) as the fraction of uninfected 
walkers. The fraction of uninfected walkers is, in fact, 
isomorphic to the survival probability of A particles in 
the reaction A + B — ► B, with equal mobility for both 
species pq,E3,p3,E6| . At late times, one expects the 



asymptotics of this model to be identical to those of the 
A + B — > diffusion- limited process with pa(0) <C Pb(0)- 
For this process it was proved rigorously by Bramson and 
Lebowitz [E3 that 



p A (t) ~ exp(-At 1/2 ) 



(2) 



in one-dimension, with unknown constant A. Therefore 
we expect identical asymptotic behavior for U (t) . The ki- 
netics of the two species annihilation model has received 
considerable attention over the past two decades and is 
largely understood |j 23 , 23j27j-E}]J . It is well appreciated 
that in this model, there exists an upper critical dimen- 
sion d c , below which spatial fluctuations in the initial 
distribution of the reactants play a significant role in the 
evolution of the density of the particles. This depen- 
dence on the microscopic fluctuations invalidates tradi- 
tional approaches such as the mean-field approximation 
which yields pit) ~ t . It was first shown by Tous- 
saint and Wilczek [^7| that, for equal initial densities 
(pa(0) — Pb(0)) of randomly diffusing A and B parti- 
cles, the walker density decays anomalously, according 
to pA (t) ~ yj pA(0)t^ d / i for d < d c = 4, assuming that 
the particles are distributed homogeneously at t = 0. In 
the case of this model we argued heuristically, on the 
basis of a toy model p] that, 



P(t) ~ expi-CiplDt) 1 /*} 



(3) 



where C is a constant and po — pa(0) + Pb(0). We ex- 
tend our study to include an investigation of the fraction 
of uninfected walkers in this model in one-dimension. 

The fraction of uninfected walkers, U(t), is related to 
the fraction of persistent sites, P(t), in the following man- 
ner. In evaluating P(t), one considers particles with dif- 
fusion constant D and addresses the probability that a 
given site x has never been visited by any walker up to 
time t. Equivalently, one could consider the site x to con- 
tain a particle with diffusion constant D' — and address 
the probability that the particle remains uninfected (i.e. 
unvisited). However, if the particle initially at position 
x is given a diffusion constant equal to that of the other 
particles in the system i.e. D' = D, then the probability 
that it "persists" up to time t is the probability that it 
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remains uninfected. This suggests that the number of un- 
infected particles is a function of the diffusion constants 
in the problem, with the limiting case D' = corre- 
sponding to the standard persistence problem. Clearly, 
in order to avoid contamination, remaining stationary is 
a more effective strategy than diffusing. Therefore, we 
anticipate that the fraction of uninfected walkers decays 
more rapidly than the fraction of uninfected (unvisited, 
i.e. persistent) sites. This faster decay of U{t) is clearly 
observed in the simulations. 

Given the similar decay of P(t) in the A% < — > %A and 
A+B — > processes, these models make a sensible combi- 
nation to study together. We discover that in both cases, 
U(t) ~ [P(t)}' t where 1 < 7 < 1.39. In fact the relation 
U(t) = k[P{t)X<, with k close to unity and 7 = 1.39(1), 
describes the data rather accurately in both models over 
the whole range of times studied, although there are in- 
dications in both models that the asymptotic value of 7 
is smaller than 1.39. We discuss the A% < — > %A model in 
Sec. II and follow with the A+B — > reaction in Sec. III. 
We conclude in Sec. IV with a discussion and summary 
of the results. 



II. THE A0 



%A MODEL 



We consider a system of non-interacting, unbiased, dif- 
fusing walkers, with random initial conditions. Our nu- 
merics are performed on a ID lattice of size L = 10 7 
with lattice constant a = 1 and periodic boundaries. At 
t = 0, N random walkers are randomly distributed on the 
lattice with a maximum of one particle per site so that 
[7(0) = 1 and P(0) = 1 - p , where p = p(0) = N/L. 
Clearly, P(0) ^ 1 is a consequence of our lattice descrip- 
tion and in the continuum limit p$ — > 0, where the mean 
distance between walkers is much larger than the lat- 
tice spacing, P(0) — > 1. Therefore, in order to place the 
fraction of persistent sites and the fraction of uninfected 
walkers on an equal footing in our numerical study, we 
define, 



Pit) 



P{t) 
P(0) 



(4) 



so that both ?7(0) and p(0) are normalized to 1. We now 
use p{t) in the description of our numerical results. 

Our model is updated using the direct method jj2| i.e. 
at each computational step, a particle is picked at random 
and shifted with probability D = 1/2 to a neighboring 
site, where D is the diffusion constant. For t > 0, any site 
which is visited by a walker becomes nonpersistent for all 
remaining time t. Mutual infection occurs when two or 
more particles simultaneously occupy a single lattice site. 
For each jump made by a particle, time t is increased by 
dt = 1/N, where N is the current number of particles 
in the system (which is constant in the present model). 
Our results are averaged over 100 independent runs. 



Dimensional analysis demands that, on a continuum, 
the fraction of unvisited sites is a function of the dimen- 
sionless combination of parameters in the problem, i.e. 
Pit) = fiplDt). Similarly, for the fract ion of uninfected 
walkers, U(t) = gip^Dt), where g is some other function. 
However, given that our numerical simulations take place 
on a discrete lattice, we expect the scaling of P(i) and 
Uit) as a function of p\P>t to become strictly true only 
in the continuum limit, po — > 0. In order to approximate 
this limit as closely as possible, we choose as our densities 
po = 0.007, 0.005 and 0.003, so that p < 1. 

The dimensionless combination of parameters, p\P>t, 
also defines a natural characteristic time scale in the 
problem, namely t* = l/p^D. Therefore one expects 
to observe the asymptotic behavior of the system only 
after a time t > i* has been reached. However, as far 
as the persistent properties are concerned, Eq. (|l|) is an 
exact relation in the continuum limit for all t. There- 
fore, we expect to observe Pit) approaching asymptopia 
very quickly in our numerical data. Clearly, in the dis- 
crete description of the model used in our simulations, 
Eq. (|l|) breaks down at the earliest times since, for very 
small t, lattice effects play a dominant role. Eq. (||), how- 
ever, which holds for U it) is a truly asymptotic relation. 
Unfortunately, the extremely fast decay of Uit) leads to 
small-number effects at very early times in our simula- 
tions and we find the limit 4>i* unattainable with good 
statistics. Therefore, we limit ourselves to studying the 
regime up to t max ~ t*/2. 

The fraction of unvisited sites decays, in the continuum 
limit, exactly as the stretched exponential form, Eq. (1) 
H . In order to make a direct comparison between the 
fraction of persistent sites, pit), and the fraction of unin- 
fected walkers, Uit), we present these two quantities on 
a log-linear plot in Fig. [I] , using ip^Dt) 1 / 2 as abscissa. 
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FIG. 1. Log-linear plot of the fraction of unvisited sites 
pit), and the fraction of uninfected walkers Uif), as a func- 
tion of ip 2 Dt) 1/2 , for the A0 < ► model. 
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Fig. |] shows clearly that the data scale as a function of 
the dimensionless quantity p^Dt. U(t) exhibits behavior 
characteristic of Eq. (||), with a faster decay rate than 
for the fraction of persistent sites i.e. A > [l&p^D /it) 1 / 2 . 
The data seem consistent with 



behavior of pit) is well described, for all t, by an exact 
lattice calculation [E5| which predicts 



U(t) ~ K\p(t)F , 



(5) 



with K ~ 1 and 7 > 1. The log-log plot of U(t) vs. p(t) 
in Fig. H indeed shows an almost linear relationship be- 
tween \nU{t) and lnp(t). In order to evaluate 7 and K 
we performed a linear regression on the data in Fig. ^ 
in the region indicated by the arrows, thereby avoiding 
initial transients associated with our lattice description, 
and the onset of statistical fluctuations between different 
runs of the data. We summarize our results for 7 and K 
in Table 1. 




-1.65 -1.15 -0.65 -0.15 

In <p(t)> 

FIG. 2. Log-log plot of the fraction of uninfected walk- 
ers, U(t), against the fraction of unvisited sites, p(t), in the 
A® < — t %A model. 



Po 


7 


K 


0.007 
0.005 
0.003 


1.391(5) 
1.391(5) 
1.391(5) 


0.990(3) 
0.991(3) 
0.991(3) 



TABLE 1. Numerical values of 7 and K, where U(t) ~ 
K\p(t)] J , for various po, in the A% < > %A model. 

The data are sufficiently good to merit closer in- 
spection. To check how far our numerics probe the 
asymptotic regime, we plot — (pyDt) -1 / 2 \np(t) and 
-(p^L»t)- 1 / 2 lnC/(i) against ]n(f%Dt) in Fig. |. The na- 
ture of this plot greatly expands the early time regime. 
Whereas the data-collapse for different densities is excel- 
lent in Fig. [l], a systematic splitting of the curves as a 
function of po is now clearly evident. This feature is a 
consequence of the discreteness of the lattice. In fact the 



p{t) = cxp[~p af(Dt/a 2 )} 



(6) 



where a is the lattice spacing. The function f{x) has 
the limiting behavior f(x) — > 2x for x — > and f{x) — > 
{l&x/ir) 1 / 2 for x — » 00, the latter reproducing the con- 
tinuum limit result, Eq. (Q). Indeed, the structure of Eq. 
(^) shows that the continuum limit, a — > at fixed t, and 
the asymptotic limit, t — > 00 at fixed a are the same. Eq. 
(||) also shows that the quantity — {p^Dt)^ 1 / 2 lnp(t) is, 
for fixed a and D, a function of t only, i.e. independent 
of po. For small t it behaves as t 1 / 2 , accounting for the 
small-i behavior in Fig. (|^). Plotting against ln(poZ3i) 
displaces the curves for different densities in the horizon- 
tal direction, as is clear in the figure. The whole region 
where the curves are split in Fig. corresponds to the top 
left-hand corner of Fig. |l|. If the data are plotted against 
hit in Fig. (||), instead of ln(pg£>i), the data for different 
densities recollapse, but the onset of the continuum limit 
at late times not as clear. 
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FIG. 3. -{plDt)- 1/2 \np(t) and -{p 2 Dt)~ 1/2 In U(t) plot- 
ted against \n{plDt) for p = 0.007, 0.005, 0.003. 

For large t, where the curves for the different densi- 
ties converge in these scaled variables, the persistence 
data clearly approach the limiting value 4/y / 7r, consistent 
with the continuum limit displayed in Eq. (|l|). By con- 
trast, Fig. H shows that all of our data for U(t) is merely 
preasymptotic. The data reach a maximum and then ap- 
pear to monotonically decrease in the region where the 
different curves converge (the continuum limit). Our re- 
sults are consistent with the recent simulations of Mehra 
and Grassberger |2(| who probe U(t) to much smaller 
values than achieved here. The slow approach to asymp- 
topia for U (t) is not understood and, although our nu- 
merics seem roughly consistent with Eq. (^) , a numerical 
determination of A appears to be a hard problem. In the 
light of Fig. ||, Fig. || and its associated results (Table 1) 
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are somewhat misleading. Clearly, over the times cov- 
ered by our simulations, 7 can not be a constant since A 
is not. 

In Fig. H we plot the effective exponent j(t) = 
ln?7(i)/lnp(£) directly as a function of time. The data 
show that 7 decreases monotonically with time even in 
the "continuum regime" where the curves for different 
densities converge. For t — > 00, 7 must tend to a constant 
if Eq. (|^) is correct. We note that since the fraction of 
uninfected walkers decreases more rapidly than the frac- 
tion of persistent sites, 7 > 1. Therefore, if the trend in 
the data continues, we can bound 7 by 1 < 7 < 1.39: 
Asymptotically, U(t) ~ [P(£)] 7 for some fixed 7 where 
1<7<1.39. 



III. THE A + B 



MODEL 
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FIG. 4. Plot of 7 vs. ln(po-Di), where 7 = In U{t)/ lnp(t), 

for the A% < > %A model. Inset: blow-up of the right-hand 

part of the plot. 

Let us consider once again the probability of being un- 
infected of a particle at site x, with diffusion constant 
D' , and diffusing particles to its left and right with dif- 
fusion constant D. Then, we suggest that the fraction of 
uninfected walkers can be expressed as, 



U(t) ~ eM-(p 2 Dt) 1/2 F(D'/D)} 



(7) 



for large t, consistent with the result of Bramson and 
Lebowitz (Eq. (||)). F(x) is some unknown function 
of the ratio of the diffusion constants in the system. 
For D' = 0, the problem reduces to the familiar study 
of unvisited sites. Therefore F(0) = (16/n) 1 / 2 and 
U(t) = P(t). The case D' = D corresponds to equally 
mobile particles, and the particle at site x becomes in- 
fected more quickly than if it had remained stationary. In 
this case F(l) = 7(16/71-) x / 2 giving U(t) ~ [P(i)] 7 where, 
according to our numerics, 1 < 7 < 1.39. Although we 
have explored the values D' = and D' = D only, we 
conjecture that Eq. ([?]) holds for any ratio of the diffusion 
constants. 



We now address the fraction of uninfected walkers, 
U(t), in the two species annihilation model. Our simula- 
tions are performed on a ID lattice of size L = 10 7 with 
periodic boundary conditions. At t — 0, exactly equal 
numbers, Na{0) = Nb(0), of A and B particles are ran- 
domly distributed on the lattice with a maximum of one 
particle per site, such that at large scales, both densities 
Pa(0) and ps(0) are initially homogeneous. We define 
Pa(0) = Pb(0) = N A (0)/L and Po = p A (0) + p B (0). 
Both species are also given the same diffusion constant 
Da — Db — D — 1/2. Our model is then updated in the 
same manner as for the A% < — ► %A model described in 
Sec. II, i.e. dt = 1/N(t), where N(t) = N A (t) + N B (t), 
is the total current number of particles in the system. 
Infection occurs when two or more particles occupy a 
single lattice site, but we impose an instantaneous reac- 
tion (A + B — * 0) so that each lattice contains only one 
type of particle. 

It is well known that in the two species annihilation 
model in one dimension, there is a effective repulsion be- 
tween the A and B particles that favors segregation into 
single species domains P, p7[pl| . This domain coarsen- 
ing leads to early time transients for both P(t) and U(t). 
Although using a high initial density serves to acceler- 
ate the systems progression into the asymptotic regime, 
in doing so, not only is the continuum limit obscured, 
but we also find that the extremely fast decay of both 
P(t) and U(t) leads to small-number effects at very early 
times. Therefore we choose po <C I . For consistency, we 
select our parameters to be the same as those studied in 
the A<D < — > process, i.e. p = 0.007, 0.005 and 0.003, 
and p\P>t max ~ 0.5 where t max is the maximum num- 
ber of time steps in the simulations. All the results are 
averaged over 100 runs. 

The format of our analysis is much the same as in the 
preceding section. We observe that an algebraic rela- 
tion, U(t) ~ [p(i)] 5 , between U(t) and p(t) also holds ap- 
proximately in this case over the range of times studied. 
Surprisingly, S seems to have a similar value to 7 in the 
A$ < — > %A model, i.e. 1 < 8 < 1.39. For the two-species 
annihilation process we suggested in our earlier paper, on 
the basis of a toy model, that the fraction of unvisited 
sites decays asymptotically according to the stretched ex- 
ponential form, Eq. (||) M . In the manner of Fig. |f|, to 
make a direct comparison between the fraction of per- 
sistent sites P(i), and the fraction of uninfected walkers 
U(t), we present these two quantities on a log- linear plot 
in Fig. ||, this time using (pgDt) 1 / 4 as abscissa, as sug- 
gested by Eq. (||). Fig. || shows clearly that the data 
scale excellently as a function of the dimensionless quan- 
tity p\B>t. We appreciate, however, that the stretched 
exponential behavior of Eq. (|J) is not particularly well 
realized. This may be due to the prevalence of large ini- 
tial transients (we presented more convincing evidence 



4 



of the exponential behavior in our earlier paper In 
particular, the asymptotic behavior represented in Eq. 
(3) is only expected for large values of PgDt, unlike Eq. 
(1) for P(t), which we have shown holds for all t on the 
continuum Nonetheless, even in the regime we have 
studied it is clear from Fig. || that the fraction of unin- 
fected walkers has a very similar form of decay to the 
fraction of unvisited sites, but with a faster decay rate. 
In Fig. ||, therefore, we plot lnU(t) against \np(t) to see 
whether there is a simple relationship between U(t) and 
p(t). 

i 
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FIG. 5. Log-linear plot of the fraction of unvisited sites 
p(t), and the fraction of uninfected walkers U(t), as a func- 
tion of p^Dt, for the A + B — > model. 
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FIG. 6. Log-log plot of the fraction of uninfected walk- 
ers, U(t), against the fraction of unvisited sites, p(t), in the 
A + B -» model. 

The approximate linearity of the data in Fig. ^ suggests 



U(t) ~ k[p(t)] s . To determine S and k, we performed a 
linear regression on the data in the region indicated by 
the arrows, thereby avoiding initial transients associated 
with our lattice description, and the onset of statistical 
fluctuations between different runs of the data. We sum- 
marize our results for S and k in Table 2. Notice that S 
and k agree remarkably well with 7 and K recorded in 
Table 1 for the A% < — > %A process. 



Po 


5 


k 


0.007 
0.005 
0.003 


1.390(5) 
1.390(5) 
1.386(5) 


0.990(3) 
0.993(3) 
0.991(3) 



TABLE 2. Numerical values of S and k, where U(t) ~ 
k[p(t)] s , for various po, in the A + B — > model. 

In analogy to Fig. 0, we plot (p§-Dt) -1//4 lnp(f) and 
(plDt)- 1 /* In £/(*) against \n(p%Dt) in Fig. @. 
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FIG. 7. -(poDt)- 1/4 \np(t) and -(plDty 1/4 In U(t) plot- 
ted against In(plDt) for p = 0.007, 0.005, 0.003. 

Early-time lattice effects, where the different densities do 
not overlap, are evident in Fig. [7]. It is possible that all 
of our data for pit) and U(t) are within the preasymp- 
totic regime, since the p(t) curve has not yet saturated to 
the late-time asymptotic values impied by Eq. (||). Note, 
however, that if the curves in Fig. |?] were asympotically 
linearly increasing, rather than eventually saturating to 
a constant, the corresponding forms for p(t) and U(t) 
would be exp [-CiplDt^l^lniplDt)], with C a constant 
(different for p and U). This form differs only by the 
logarithmic term from that predicted in [|| on the ba- 
sis of a toy model. It should be noted that neither our 
present data nor that presented in Q can definitively 
rule out such logarithmic corrections, and the asympotic 
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form of p(t), as well as U(t), cannot be regarded as be- 
ing definitively established. Throughout the preasymp- 
totic regime, however, the relation U(t) = K[p(t)] s , with 
S = 1.39(1), holds rather well. In Fig. [| we plot the effec- 
tive exponent 8(t) = In U(t) / \np(t) directly as a function 
of \n{plDt). 




-12.75 -7.75 -2.75 

In (p/Dt) 

FIG. 8. Plot of 8 vs. \nt, where S = In U(t)/ lnp(t), for the 
A + B — > model. Inset: blow-up of the right-hand part of 
the data. 

Fig. H is remarkably similar to Fig. [| for the A% < — > A 
process. Indeed, except for the earliest times, the cor- 
responding curves differ by less by 1% everywhere. In a 
similar way we can bound the asymptotic S (assuming 
that the late-time trend evident in Fig. || continues) by 
1 < 5 < 1.39. Since we have argued elsewhere ||, using 
a toy model, that the asymptotic behavior of P(t) is de- 
scribed by Eq. (j|), we conjecture that the asymptotics of 
U{t) are given by 



U(t) ~ expi-SCiplDt) 1 / 4 } 



(8) 



where C is a constant and 1 < 8 < 1.39. For the 
A% < — ► %A model, we expressed [Eq. (^)] U(t) in terms 
of the ratio D' /D, where D' is the diffusion constant of 
a tagged particle, D is the diffusion constant of the other 
particles, and U(t) is the probability that the tagged par- 
ticle remains uninfected at time t. The analogous form 
for the A + B — ► process is naturally, 



U(t)~eM~(p 2 oDt) 1/4 G(D'/D)] 



(9) 



where G(x) is some unknown function. As before, the 
case D' = corresponds to the problem of unvisited sites, 
so that G(0) = C and U(t) = P(t). For D' = D, the case 
studied above, our conjecture implies G(l) = SC. How- 
ever, we stress that a logarithmic correction of the form 
discussed above cannot be ruled out by the data. The 
main motivation for Eq. ([)]) is the toy model presented 
in KJ], applied to the case D' = 0. 



IV. DISCUSSION AND SUMMARY 



In this paper we have studied two distinct models: (i) a 
one-dimensional system of non-interacting, unbiased dif- 
fusing walkers A% < — ► 0^4, with random initial condi- 
tions and (ii) the A + B — * diffusion limited process 
in one dimension with equal initial densities of A and B 
particles distributed homogeneously at t = 0. We have 
reduced the familiar study of persistent sites P(t), to 
a limiting case in the study of uninfected walkers and 
shown numerically that in both models, 



U(t) ~ [P(t)]T , 



(10) 



where 7 ~ 1.39 for the time regimes covered by our data 
and 1 < 7 < 1.39 for asymptotically large time. The 
common feature of the A% < — ► %A and A + B — > 
processes is that the fraction of unvisited sites decays 
according to a stretched exponential (possibly with loga- 
rithmic corrections for A + B — > 0). In our earlier paper 
||, we argued heuristically, on the basis of a toy model, 
that if the density of walkers decays asymptotically as 
t~ a , then the fraction of unvisited sites P(t) decays with 
a stretched exponential form if a < 1/2. We showed that 



P(t) - exp-[A(a)plDt] 



(1/2)- 



a<l/2 (11) 



where A is some function of a. This result is consistent 
with the results for models studied here. Remarkably, 
Eq. @ seems to hold for both the A0 < — ► 0A and 
the A + B — > processes, with approximately the same 
bounds on 7. Furthermore, the effective time-dependent 
value of 7 is very similar for both models (see Figs. ^ and 
||). An obvious goal for the future is to try to obtain a 
theoretical understanding of this simple relationship, and 
test it against other models with a < 1/2. 

For both models studied here we expressed the prob- 
ability that a particle at site x is uninfected in terms 
of its diffusion constant D' and the diffusion constant 
D of the other particles [Eqs. (j7p and (||)], noting that 
P(t) corresponds to the limit D' = 0. Using Eq. ( pd| ) , 
we can generalize this description to any system of ran- 
domly diffusing particles for which the density decays as 
~ t~ a with a < 1/2. We conjecture that, in this case, the 
fraction of uninfected walkers can be expressed, asymp- 
totically, as 



U(t) ~ exp[-(pgDt) 



(1/2)- 



l R(D'/D,a)} 



(12) 



with U(D' = 0) = P(t) and R some unknown func- 
tion. Clearly, there is much scope to test the validity 
of this form since we have only considered a = 0, 1/4 
and D' = 0, D. Note that for systems with a > 1/2, our 
toy model predicts [|| that P(t) approaches a non-zero 
constant, while U(t) certainly behaves differently since 
it is bounded above by p(t) and must therefore vanish 



G 



for large t. More interesting, however, is the borderline 
case where the particle density falls off like i -1 / 2 . In 
this case, our toy model predicts that P(t) |(| decays as 
a power-law, P(t) ~ t~ e . A system where such behav- 
ior is observed is the well studied g-state Potts model 
|,|jl|Jl|Jl|. The results of our study of uninfected 
walkers in this model follows in a separate paper j33[|. 
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